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Recently, we developed a mean-field-type framework which treats the correlation induced by the 

l> . 

tensor force. To exploit the tensor correlation we introduce single-particle states with the parity 

\Q ■ 

and charge mixing. To make a total wave function have a definite charge number and a good parity 



the charge number and parity projections are performed. Taking a variation of the projected wave 
function with respect to single-particle states a Hartree-Fock-like equation, the charge- and parity- 
projected Hartree-Fock equation, is obtained. In the charge- and parity-projected Hartree-Fock 
method, we solve the equation selfconsistently. In this paper we extend the charge- and parity- 
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projected Hartree-Fock method to include a three-body force, which is important to reproduce the 
saturation property of nuclei in mean-field frameworks. We apply the charge- and parity-projected 
Hartree-Fock method to sub-closed-shell oxygen isotopes ( 14 0, 16 0, 22 0, 24 0, and 28 0) to study 
the effect of the tenor correlation and its dependence on neutron numbers. We obtain reasonable 
binding energies and matter radii for these nuclei. It is found that relatively large energy gains 
come from the tensor force in these isotopes and there is the blocking effect by occupied neutron 
orbits on the tensor correlation. 
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I. INTRODUCTION 



The tensor force plays important roles in nuclear structure. The study of nuclear matter 
using the Brueckner theory showed that the tensor force is closely related to the binding 
mechanism and the saturation property of nuclear matter [l| . The almost exact calculations 
with very large model spaces exhibit a large attractive energy comes from the tensor force 
. The tensor force is inferred_to_be responsible for about a half of the single-particle 



spin-orbit splitting in light nuclei 



4. 



Recently due to the development of experimental techniques, we can access to various 
kinds of unstable nuclei experimentally. Those experiments have been revealing that the 
shell structures of unstable nuclei may change from those of stable nuclei [tUs!]- Considering 
the importance of the tensor force in nuclear structure, the tensor force probably has an 
effect on such structure changes of nuclei jjj. Therefore, the study of the effect of the tensor 
force in neutron-rich nuclei is interesting and important. 

To study the effect of the tensor force in a relatively large mass region in the nuclear 
chart includingunstable nuclei, we have developed a framework based on a mean-field-type 
model One of the most important tensor correlations in closed-shell nuclei is 

a 2-particle-2-hole (2p-2h) correlation. In a usual Hartree-Fock calculation, the correlation 
induced by the tensor force can be treated restrictively and the 2p-2h correlation is hard to 
be handled. The effect of the tensor force is thought to be included in other kinds of forces 
like the central, LS, and density-dependent forces in the usual Hartree-Fock calculations. To 
treat the tensor force directly, we introduce a single-particle state with the parity and charge 
mixing considering the pseudoscalar and isovector characters of the pion, which mediates 
the tensor force [10|, [ll|. Because a total wave function made from such single-particle 
states with the parity and charge mixing does not have a good parity and a definite charge 
number, the parity and charge number projections are performed before variation 
We call this method the charge- and parity-projected Hartree-Fock (CPPHF) method. In 
the previous studies we applied the CPPHF method to the alpha particle and showed that 
the tensor correlation can be exploited in the CPPHF method. The CPPHF method was 
also applied to 8 Be to study the effect of the tensor force on alpha clustering 

There are other attempts to treat the tensor correlation by expanding usual model spaces 
like a mean- field model Q], a shell model jfj, the antisymmetrized molecular dynamics 
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(AMD) [16]. Those studies including ours showed the importance of the 2p-2h configuration 
mixing and high-momentum components in single-particle states, which are not treated in 
usual model space calculations. Neff and his collaborators took a different approach to treat 
the tensor correlation They used the unitary correlation operator method (UCOM) 

to make an effective interaction in a moderate model space. In their effective interaction, 
the attractive correlation by the tensor force is included in other forces like the central and 
LS forces. They performed the calculations up to the second order perturbation based on 
Hartree-Fock calculations using their effective interaction and obtained a nice agreement of 

n q 

binding energies over the whole mass region 18]. Otsuka and his collaborators [9j] showed 
that a particle-hole (p-h) correlation by the tensor force is also important and changes single- 
particle spin-orbit splittings of neutron (proton) orbits with proton (neutron) numbers. The 
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effect was inferred in old Hartree-Fock calculations 

In the present paper, we apply the CPPHF method to sub-closed-shell oxygen isotopes, 
1 O, 16 0, 22 0, 24 0, and 28 0, which are assumed to have sub-closed-shell structures for 
neutron orbits up to 0^3/2, Opi/2, lsi/2, and 0d 3 /2 respectively, to see the dependence 

of the correlation induced by the tensor force on neutron numbers. We extend the CPPHF 
method to treat a three-body force, which is needed to reproduce the saturation property 
of nuclei with relatively large mass numbers. In Section |H] we explain the CPPHF method 
with a three-body force. In Section ITTT1 the results of the CPPHF method are presented. In 
Section HH we summarize the paper. 



II. CHARGE- AND PARITY-PROJECTED HARTREE-FOCK METHOD WITH 
A THREE-BODY FORCE 

In this section we formulate the charge- and parity-projected Hartree-Fock (CPPHF) 
method in the case where a three-body force exists. A Hamiltonian for an A-body system 
with the two-body and three-body forces can be written as 

A A A 

H = ^t(x a ) + v {2) {x a ,x b )+ ^ 3 \x a ,x b ,x c ), (1) 

a=l a>b=l a>b>c=l 

where t, i)^, and are one-body, two-body, and three-body operators respectively, x's 
are coordinates including spin and isospin. In the CPPHF method, we assume as single- 
particle states the ones with the parity and charge mixing. It means each single-particle 
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wave function has both positive-parity and negative-parity components and both proton 
and neutron components. The single-particle wave function which consists of these four 
components has the following form, 

M*)= E E VW*»- (2) 

p a =± t xa =±l/2 

In the above equation, p a denotes parities, + for positive parity and — for negative parity, 
and t za denotes isospins, +1/2 for proton and —1/2 for neutron. In the CPPHF method, 
we take as a wave function of an A-body system a Slater determinant which consists of the 
single-particle states with the parity and charge mixing, 

1 A 

Here, A is the antisymmetrization operator. Because ^/ mtr does not have a good parity and 
a definite charge number, we need to perform the projection operators of parity (±) and 
charge number (Z) on \l/ mtr to obtain the wave function with a good parity and a definite 
charge number; 

y(±;Z) = fip(±}fic(z)m™ tT . (4) 

Here, V p (±) is the parity-projection operator, where V p (+) projects out the positive parity 
state and P p (— ) projects out the negative parity one. V C (Z) is the charge-number-projection 
operator, which projects out the wave function with a charge number Z. Therefore, \I>( ±;Z ) 
has a good parity (±) and a definite charge number (Z). The parity projection operator 
P p (±) is defined as 

P p (±) = ^ (P = f[p^> (5) 

where the total parity operator P is the product of the parity operator p a for each single- 
particle state. The charge projection operator V C (Z) is defined as 

(6) 

where Z is the charge number operator, which is the sum of the single-particle proton 
projection operator (1 + r^)/2, and the charge-rotation operator is defined as C{9) = e %ze . 




We take the expectation value for the Hamiltonian H with the projected wave function 
and obtain the energy functional, 

E (±-,z) _(y {± ' Z) \H\y {± '' Z) ) _ (ty int *\H\V p (±)V c (Z)^' mtv ) 

(V[/(±;^)|V|>(±;Z)) /^intr|pp^pc^Wntr\ 



± d9e- iZd (n(°) (9) ± n( p ) (0)) ' 

The denominator in the right-hand side of the above equation is the normalization of the 
total wave function, 

n (± ; z) = ^(±;z)|^(± ; z)) = _L r d ee~ iZ6 (n<°>(0) ±^(0)) . (8) 

*/ 

Here, n^(9) is the determinant of the norm matrix between the original single-particle wave 
functions ip aa and the charge- rotated single-particle wave functions ip aa (9). n^(9) is the 
determinant of the norm matrix between the original single-particle wave functions ip aa and 
the parity-inverted and charge- rotated single-particle wave functions ip^(9). 

n^{9) = (* intr |C(0)|* intr > = det B^(9) (B^(9) ab = (VU^M*))), 

(9) = (¥^\PC(9W ntI ) = detB( p \9) (B^(9) ab = (VU^S? (*)>)• (9) 



(7) 



n 



(P) 



The charge-rotated wave function ip aa (x b ; 9) and the parity-inverted and charge-rotated wave 
function ipa](x b ] @) are defined as 

=e ie i) Pa=+)tza=l/2 {x b ) + e ie ip Pa= „ itza=1/2 {x b ) 

+ *Pp a =+,t za =-l/2(Xb) + 1pp a =-,t za =-l/2(Xb), (10) 

=e t8 ij pa=+jtza=1/ 2(x b ) - e t6 ij pa= _ ttza=1/2 (x b ) 

+ ^ Pa =+,t za =~l/2(x b ) - i>p a =-,t za =-l/2(x b ), (11) 

where p b is the single-particle parity operator in (jSJ) and (1 + r b )/2 is the single-particle 
proton projection operator in (JHJ). 

The numerator in the right-hand side of ((Zj) is the unnormalized total energy, 

(*(±i*)|jf|*(±i3)) = _L r d9e- lZd (E<-°\9)±E<P\e)) . (12) 

J 
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E(°>(9) in the right-hand side of (J12)) has a similar form as a simple Hartree-Fock energy but 
the single-particle wave functions in the ket are modified by the charge rotation, 

E (0) (9) = (¥ ntI \HC(9)\¥ ntI ) 

A A 

a=l a>b=l 

A 

+ (^Ma^ltaM^aM^M)- (13) 



a>b>c=l 



Here, ip aa (x] 0) is the superposition of 4> aa (x] 0) weighted by the inverse of the charge-rotated 



norm matrix (B^(9) 1 



ba ■ 



A 



KM) = E^( X ^)( 5(0) (^V ( 14 ) 

6=1 

This summation for ip ab (x; 9) comes from the antisymmetrization of the total wave function. 
The hats in the last two terms represent the antisymmetrization as 

+ 4> ac (9)i> aa (9)i> ab (9) - $ aa (9)^ ac (9)i> ab (9) 

E^°\9 = 0) reduces to a simple Hartree-Fock energy. £ (p) (6>) in the right-hand side of JEJ) 
has a similar form as E^ \9) but tjj aa (9) , s are replaced by ipa](9) , s, 

£ (p) (0) = (¥ ntI \HPC{9)\¥ ntv ) 

A A 

=E^«-W£ ) W>+ E (1>cJ>a b \*\$$!(e)$%(0)) 

a=l a >b=l 

A 

+ E (^JaM^ { £(0)^ b \0)^(9)). (17) 



a>fe>c=l 

Here, ipa](x;9) is the sum of il)a](x;9) weighted by the inverse of the parity-inverted and 
charge-rotated norm matrix (B^(9)~ 1 )i Ja , 

A 

= E^^^c^Cf)- 1 )^ (is) 

6=1 



We then take the variation of E^ ±]Z ^ with respect to a single-particle wave function ip aa , 

A 



E (±;z) _ e 6c (^ 6 |Va c > \ = 0. 

b,c=l 



(19) 



The Lagrange multiplier e a b is introduced to guarantee the ortho-normalization of single- 
particle wave functions, (Vwl^at) = ^a a ,a b - As the result, we obtain the following Hartree- 
Fock-like equation with the charge and parity projections (the CPPHF equation) for each 



47T 



/•27T 




/ d9e~ ize 


n<°>(0){ 


Jo 



6=1 



±n( p ) (5) { (/> (1)(p) (*«; 0) + /> (2)(p) (x a ; 0) + ^ (x a ; 0)) 

- - EW(0))jV>(x a] 9) - £r£\0)jV>(x a] 9) 

b=i 

A 

= n {± ' z) ^2e ab ip ab (x a ), 
b=i 

where a = 1, 2, . . . , A. Here, r]^(9) and i]^(9) are defined as follows, 

A ______ 

V l 2W) ^aj\i; ab (9)) +J2(^^v {2) \^a b (0)^M) 



c=l 



c,d=l 



c=l 



(20) 



(21) 



(22) 



c,d=l 



M 1 ), /i^ 2 \ and ft/ 3 ) are the single-particle operators originated from the one-body, two-body, 



and three-body operators, which are defined as following, 



A 
6=1 

- (^v {2 \x a )\^ a M)i^ b (x a -9) 
h (3) ^ aa (x a ;9) =- ^[{^ ah i) a ^\x a )\^ ab {9)^ a X9))i^ aa {^ 



(23) 



(24) 



6,c=l 



+ (V'a^cl^ (3) (^)|4 (^)^ a W}l,2V5 Ob l^ai "J 

+ (<^J^ (3 H^)|^ WX b W)i,2^ c (x a ^)}. (25) 

The expressions of h^^ p \ h^ 2 ^ p \ and are obtained from those of h^ 2 \ and Z?/ 3 -' by 

replacing ^ a {x, 9) with 9). The notations for the integration of the two-body matrix 

elements, 



(ip ab \v {2) {x a )\i! ac ) 1 = / dx 1 i/)l lb (x 1 )v^\x a ,Xi)'4) ac (x 1 ) 



(26) 



and for that of the three-body matrix elements, 

(lpa b 1pa c I V (3) (X a ) | ^a d 1pa e ) 1 , 2 

cfei y rfx 2 V'a 6 (^i)V'iJ^2)w (3) ( a; a,a;i,a;2)^ ad (a;i)V' ae (a;2) (27) 

are introduced. The system of the coupled equations (|2T?j) for a = 1, • • • , A are solved 
selfconsistently. 

We give here the expressions for the expectation value of the kinetic energy (T)^'^ with 
the center of mass correction, that of the two-body potential energy (v 

(2)^(± ; z) ; and that of 

three-body potential energy (■0( 3 ))( ± ' z ) for the later convenience. 

1 /-27T 



4 7m (±;Z) 



d9e 



-ize 



n 



(0) 



h 2 



. 1=1 

.4 



a>b=l 
A 



h 2 



a=l 



a>fe=l 



AM 



(28) 
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r-27T A 



a>b=l 

±n( p )(^)(^^J^ (2) l^ P a ) (^ ) W) \, (29) 
(^ (3) ) (±;Z) = i ^y r^e-^ ^ (^ (0 H^)(^>^^J^ (3 ¥« a WCw^ c W) 

71 ^° a>6>c=l I 

± (9) (^>a^a c |^ (3) |$£ (e)^)$<£ (9)) J . (30) 

III. APPLICATIONS TO SUB-CLOSED-SHELL OXYGEN ISOTOPES 

In this section, we apply the charge- and parity-projected Hartree-Fock (CPPHF) method 
formulated in the last section to the ground states of the sub-closed oxygen isotopes, 14 0, 
16 0, 22 0, 24 and 28 0. We assume the parities of these nuclei in the ground states are 
positive. Those nuclei are assumed to have the closed shell up to 0pi/ 2 for proton and 
sub-closed or closed shells up to 0p 3 / 2 ( 14 0), 0pi/ 2 ( 16 0), 0d 5 / 2 ( 22 0), ls 1 / 2 ( 24 0) and 0d 3 / 2 
( 28 0) for neutron. Although 28 is know to be unbound from the experiment, we calculate 
the nucleus to study the shell-configuration dependence of the contribution from the tensor 
force theoretically. We assume the spherical symmetry. In this case, only the total angular 
momentum j is a good quantum number of a single-particle state, because parities and 
charges are mixed in intrinsic single-particle states. An intrinsic wave function can be 
written in the following form, 

^=a n n n m 

0<J<i max -j<m<j l<nj<n™ ax 

A single-particle wave function ip nj jm is composed of four components, proton and positive 
parity, proton and negative parity, neutron and positive parity, and neutron and negative 
parity, 

il> ni im(x)= fei + t,(r)^ i+m (n)C(t,) + ^(.t^^^l^Cfe)). (32) 
t*=±i 

Here, ^(m(^) is the eigenfunction of the total angular momentum j = I + s, ((t z ) is the 
isospin wave function with t z — 1/2 for proton and t z — —1/2 for neutron. l + and l_ are the 



orbital angular momentum with positive parity and negative parity respectively. For example 
Z+ = and I- — 1 for j = 1/2, l + = 2 and /_ = 1 for j = 3/2, and so on. Eq. ()32j) indicates 
that in the present calculation assuming the spherical symmetry only the correlations among 
the same j orbits can be treated. It is a limitation of the CPPHF method with the spherical 
symmetry. For the calculation of 16 four single-particle states with j = 1/2 and two states 
with j = 3/2 are included. They correspond to vrsi/2, fSi/2, and upy 2 for j = 1/2 and 

7rp3/2 and z/p 3 / 2 for j = 3/2. Because parities and charges are mixed in single-particle states 
such a classification is approximately valid. For the calculation of 14 one state with 
is subtracted and for the calculation of 22 0, 24 and 28 new states with j = 5/2, j = 1/2 
and j = 3/2 are added one by one. 

In the present study we use the modified Volkov force No. 1 (MV1 force) |2l| for the 
central potential and the G3RS force 



23] for the non-central forces. The MV1 force is the 
modified version of the Volkov force No. 1 j^] and includes the 5-f unction- type three-body 
force, 



-,(3), 



x a , x b , x c ) = t 3 5(x a - x b )5(x b - x c ). (33) 



The Majorana parameter in the MV1 force is fixed to 0.6. The G3RS is determined from the 
nucleon-nucleon scattering data. The effect of the tensor force is effectively included in the 
MV1 force because the MV1 force is determined so as to reproduce the binding energy of 
16 in the absence of the tensor force. The effect of the tensor force is thought to appear in 
the 3 E channel of the central force as attraction. Therefore we multiply the attraction part 
in the 3 E channel of the central force by xq- We also multiply the three-body-force part by 
X3B- The (5-function-type three-body force in Eq. (|33|) reduces to the density-dependent two- 
body force, \Hp( Xa 1 2 Xb ) l+P ^ aV) 5{x a — x b ), for the wave functions of even-even nuclei with 
time-reversal symmetry in the Hartree-Fock level j^J). P a is the spin-exchange operator and 
p is a single-particle density. In this case the density-dependent force only acts on the 3 E 
channel. In the MV1 force t 3 is positive and the density-dependent force has the repulsive 
effect on the 3 E channel. The LS forces determined from the NN scattering data is usually 
weak to be used in the mean filed (Hartree-Fock) calculation. Hence, we multiply the LS 
force by 2. In this case the strength of the LS force is comparable to those adopted in the 
Skyrme forces and the Gogny forces j^ . 

We also multiply the T\ ■ r 2 part of the tensor force a numerical factor x^ as in the 
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previous study [l^, because the CPPHF method is a mean- field-type calculation and can 
only take into account the correlations induced by limited couplings among single-particle 
states. Actually, in the spherical symmetry the 2-particle-2-hole (2p-2h) correlation which 
can be treated in the CPPHF method is like (j Pl j P2 jh^jh2) with j pi = jhi and j P2 
However 2p-2h configurations with j Pl ^ jf ll and j P2 7^ jh 2 are also important 
Furthermore, other effects may enhance the tensor correlation as mentioned in our previous 
paper Q]. To take into account such effects effectively, we take (the strong tensor 

force case) in addition to xt=1-0 (the normal tensor force case), xq and X3B are determined 
to reproduce the binding energy and the charge radius of 16 for each xt- 

We expand single-particle wave functions in the Gaussian basis as in the previous study 



12j. The number of the Gaussian basis used is 10 for each orbit with the minimum range 
0.5 fm and the maximum range 10 fm. The CPPHF equation (j2*Uj) is solved by the gradient 
or the damped gradient method 2^|. The convergence of the calculation is quite slow for 
the case with a large difference between a proton number (Z) and a neutron number (N). 

n 

as in 28 0. To remedy it, the quadratic constraint potential term [2£| for Z, 

(^ intr |-(Z-Z) 2 |^ intr ), (34) 
2 

is added to the energy functional (J7J). The addition of the constrained potential makes the 
convergence faster. The value of A is taken as 1000 MeV. 



A. Results for lb O 

We first take 16 as a typical example and show the effect of the tensor force in the 
CPPHF method. In Table U] the results for 16 in the Hartree-Fock (HF) and the charge- 
and parity-projected Hartree-Fock (CPPHF) schemes are shown. In the HF scheme we fix 
Xt, xq, and x^b to 1.0. In the CPPHF scheme xq and x^b ar e 1.025 and 1.25 for the normal 
tensor force (xt = 1-0) case, and 1.040 and 1.55 for the strong tensor force (xt = 1-5) 
case. The root-mean-square charge radius R c is calculated from the proton root-mean- 
square radius R p as R c = ^/R^ + 0.64. This approximation for R c corresponds to assume 
the charge radius of proton as 0.80 fm. In the HF calculation the expectation value for the 
potential energy from the tensor force Vt is negligibly small. If we perform the charge and 
parity projection before variation (the CPPHF case), Vt comes out to be a sizable value. It 
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TABLE I: Results for 16 in the HF and CPPHF method, xt is a numerical factor multiplied 
to the T\ ■ T2 part of the tensor force. E and T are the total energy and the total kinetic energy 
respectively. Vc, V3B, Vr, Vls and Vc ou i are the potential energies form the central, the three-body 
force, the tensor force, the LS force, and the Coulomb force respectively. Those are give in the unit 
of MeV. R c is the root-mean-square charge radius in the unit of fm. 





x T 


E 


T 


V c + V m + Vcoui 


V T 


Vls 


R c 


HF 


1.0 


-124.1 


230.0 


-353.2 


0.0 


-0.9 


2.73 


CPPHF 


1.0 


-127.1 


237.1 


-351.6 


-11.7 


-1.0 


2.73 


CPPHF 


1.5 


-127.6 


253.9 


-342.2 


-38.3 


-1.0 


2.73 



becomes about 10 MeV for the normal tensor force case and about 40 MeV for the strong 
tensor force case. This result indicates that the CPPHF method is effective to treat the 
correlation from the tensor force. In the CPPHF cases the kinetic energy T becomes larger 
than in the HF case. In the CPPHF scheme, to gain the tensor correlation energy the 
opposite-parity components compared to the simple shell-model picture have to be mixed 
into single-particle states. This mixing causes an over shell correlation and, as the result, 
the kinetic energy becomes larger. A similar tendency is also observed in the alpha particle 



case 



12 



In Figs. H 121 and |3] the intrinsic single-particle wave functions in 16 with the strong 
tensor force (x^ = 1.5) are plotted. The wave functions plotted have proton components 
as dominant ones. The wave function in Fig. Qhas a proton S1/2 component as a dominant 
component. The mixing probabilities for negative parity (pi/2) and neutron are 19.5% and 
27.4% respectively. From the figure you can see that the spread of the py% components are 
smaller than the Sy2 ones. It indicates that the opposite-parity components induced by the 
tensor force have high-momentum components. The wave function which has a proton p 3 / 2 
component as a dominant one is plotted in Fig. [2] and the shrinkage of (I3/2 components 
compared to P3/2 ones are clearly seen also. In the wave function which has a proton pi/% 
component as a dominant one, which is plotted in Fig. El the shrinkage is not so clear 
compared to the previous two cases, probably because of the orthogonality condition to the 
wave function to the first j = 1/2 state in Fig. ^ 

In the alpha particle case, p\li components mixing into Sx/2 ones are also compact in size 
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FIG. 1: (Color online) Intrinsic single-particle wave function of the j = 1/2 state with an si/ 2 
proton component as a dominant one in 16 O as a function of the radial distance R. The solid, 
dashed, dotted and dashed-dotted curves correspond Sij 2 proton, p 1 / 2 proton, sj/ 2 neutron, and 
Pi/2 neutron respectively. P(— ) and P (u) are the mixing probabilities of the negative-party and 
the neutron components respectively. 

12 , 13, 3|. The importance of this shrinkage is confirmed in a s 
the antisymmetrized molecular dynamics (AMD) calculation 
infer that the shrinkage of mixing single-particle wave functions is generally important in 
heavier-mass region. There are also wave functions with a neutron component as a dominant 



lell model calculation lyl and 
la], too. The present result 



one. The general tendency just mentioned above is almost the same if proton and neutron 
are interchanged. 

In Fig.HJthe densities for 16 in the HF and CPPHF calculations are shown. Because the 
tensor force induces the opposite-parity components with narrower widths in single-particle 
wave functions, the densities are depleted in the middle in the CPPHF calculations compared 
to the one in the HF calculation. This effect is larger for the case with the strong tensor 
force as expected. To see the effect of the tensor correlation more clearly, in Fig.EJthe charge 
form factors are plotted as a function of the momentum squared. From the figure, higher- 
momentum components appear in the CPPHF calculation. It indicates that the tensor force 
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FIG. 2: (Color online) Intrinsic single-particle wave function of the j = 3/2 state with a P3/2 
proton component as a dominant one in 16 as a function of the radial distance R. The solid, 
dashed, dotted and dashed-dotted curves correspond d 3 / 2 proton, p s / 2 proton, d 3 / 2 neutron, and 
neutron respectively. P(+) and P{y) are the mixing probabilities of the positive-party and 
the neutron components respectively. 

enhances the charge form factor in a high-momentum region. The short-range correlation, 
which is not treated properly here, should have a contribution to the charge form factor in 
the high-momentum region. Hence, we need further investigation to compare the present 
result of the charge form factor in the CPPHF method with the experimental data. The 
enhancement of the charge form factor in a high-momentum region is also found in 4 He in 
the calculation with the charge- and parity-projected relativistic mean field model jl3|. 



B. Results for the oxygen isotopes 

In this subsection we show the results for the sub-closed-shell oxygen isotopes. In Fig. ® 
the results for the binding energies per particle in the Hartree-Fock (HF) calculation (the 
circle symbols), the charge- and parity-projected Hartree-Fock (CPPHF) calculation with 
the normal tensor force (xt=1.0) (the triangle symbols) and the CPPHF calculation with 
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FIG. 3: (Color online) Intrinsic single-particle wave function of the j = 1/2 state with a p\/2 
proton component as a dominant one in 16 as a function of the radial distance R. The solid, 
dashed, dotted and dashed-dotted curves correspond si/ 2 proton, pi/ 2 proton, si/ 2 neutron, and 
Pi/2 neutron respectively. P{+) and P{y) are the mixing probabilities of the positive-party and 
the neutron components respectively. 

the strong tensor force (xt=1.5) (the diamond symbols) are shown. The experimental data 



(the square symbols) 



are also plotted. The general tendency is reproduced in our result, 



although the agreement with the experimental data is not as good as the Hartree-Fock-type 



calculations 125 



31 



32 



33 



, 3J|. The HF calculation with the MV1 force underestimates the 
binding energy of 14 and 22 a little bit largely. It indicates that if we adopt our calcu- 
lation on the central and the density-dependent forces in the recent sophisticated effective 
interaction, the agreement with the experimental results should become better. There is 
an ambiguity in the treatment of the density- dependent force when we perform the parity 
and charge projections, because the density- dependent force cannot be written in a simple 
two-body operator form. It causes a difficulty when we use the density-dependent force in 
the CPPHF calculation. The optimization of the central force and the management of the 
density-dependent force will be our future problems. The agreement with the experimental 
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FIG. 4: (Color online) Densities for 16 in the Hartree-Fock (HF) and the charge- and parity- 
projected Hartree-Fock (CPPHF) methods as a function of the radial distance (R). The solid, the 
dashed, and the dotted curves correspond to the HF calculation, the CPPHF calculation with the 
normal tensor force (xt = 1-0), and the CPPHF calculation with the strong tensor force (xt = 1.5). 

data is good for 14 in the CPPHF method with the strong tensor force. In this case quite 
a large attractive potential energy comes from the tensor force as shown below. 

In Fig. the results for the root-mean-square matter radii (-R m ) are plotted with the 
experimental data [30]. We use the same symbol for each case in Fig. |U1 Except for 14 
the results for all the three cases are almost the same and reproduce the experimental data 
well. For 14 the CPPHF calculation with the strong tensor force gives a smaller matter 
radius compared to the other two calculations and the agreement with the experimental data 
becomes better in the strong tensor force case. The reduction of the radius is caused by the 
shrinkage of the opposite parity component in a single-particle wave function induced by the 
tensor correlation. Such an effect cannot be treated in simple Hartree-Fock calculations. 

To see the effect of the tensor force on the binding mechanism in the oxygen isotopes 
in Fig. |H1 the results for the total kinetic energy T (the dashed lines) and the sum of the 
potential energies from the central, the three-body and the Coulomb forces Vc + V3B + Vc ou i 
(the solid lines) are plotted. We also show in Fig. EH the result for the potential energies 
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FIG. 5: (Color online) Absolute values of the charge form factor for 16 in the Hartree-Fock 
(HF) and the charge- and parity-projected Hartree-Fock (CPPHF) methods as a function of the 
momentum squared (q 2 ). The solid, the dashed, and the dotted curves correspond to the HF 
calculation, the CPPHF calculation with the normal tensor force (xt = 1.0), and the CPPHF 
calculation with the strong tensor force (xt = 1.5). 

from the tensor force Vr (the solid lines) and the LS force Vls (the dashed line). The same 
symbols are used for the HF calculation and the CPPHF calculations with the normal and 
strong tensor forces as in Fig. All values are divided by mass numbers to make a isotope 
dependence clear. The total kinetic energy and the sum of the potential energies from the 
central, three-body and Coulomb forces show a volume-like behavior. The kinetic energy 
for 14 in the CPPHF calculation with the strong tensor force is larger than the other two 
cases. The increase of the kinetic energy is also caused by the strong tensor correlation 
in 14 0, because to gain the correlation energy from the tensor force the opposite-parity 
components must mix into single-particle states and the opposite-parity components have 
larger kinetic energy. 

The potential energies from the LS force behave in almost the same manner for the three 
cases. It becomes attractive if neutron shells are jj-closed and negligibly small if neutron 
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FIG. 6: (Color online) Binding energies per particle with minus sign for sub-closed-shell oxygen 
isotopes. The horizontal line indicates mass numbers. The circle symbols correspond to the 
Hartree-Fock (HF) calculation, the triangle ones to the charge- and parity projected Hartree-Fock 
(CPPHF) calculation with the normal tensor force (xt = 1-0), and the diamond ones to the CPPHF 
calculation with the strong tensor force (xt = 1.5). The square symbols indicate the experimental 



shells are LS-closed. The potential energy from the tensor force becomes weakly repulsive 
in the HF calculation for all the oxygen isotopes. In the CPPHF calculations the potential 
energies from the tensor force become attractive for all the oxygen isotopes. In contrast to 
the LS potential energies, the tensor potential energies have sizable values in LS-closed shell 
nuclei like 16 and 28 0. The attraction from the tensor force is the same order that from 
the LS force even in the case with the normal tensor force. In CPPHF case Vr becomes 
maximum in 14 and decreases with the mass number. For the strong tensor force case, the 
attractive energy from the tensor force becomes larger as expected. The attractive energy 
is quite large for 14 in the CPPHF calculation with the strong tensor force. For other 
isotopes the attractive energies from the tensor force are small and do not change so much 
with the mass number. 

In Fig. El the probabilities of the mixing of the opposite parity components P m i x in 



data 29]. 
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FIG. 7: (Color online) Root-mean-square matter radii for sub-closed-shell oxygen isotopes. The 
horizontal line indicates mass numbers. The circle symbols correspond to the Hartree-Fock (HF) 
calculation, the triangle ones to the charge- and parity projected Hartree-Fock (CPPHF) calculation 
with the normal tensor force (xt = 1-0), and the diamond ones to the CPPHF calculation with 
the strong tensor force (xt = 1.5). The square symbols indicate the experimental data with error 
bars 

the single-particle states with neutron components as the dominant ones in the result of 
the CPPHF calculation with the normal tensor force are shown. The result for the strong 
tensor case shows almost the same tendency. In the usual shell-model classification the 
first j—1/2 state (si/2 dominant), the first j=3/2 state (^3/2 dominant), the second j—1/2 
state (pi/2 dominant), the first j=5/2 state (d 5 / 2 dominant), the third j = 1 / 2 state (s 1 / 2 
dominant), and the second j=3/2 state ((^3/2 dominant) correspond to Osi/2, 0^3/2, Opi/2, 
0g?5/2, lsi/2, and 0d 3 /2 respectively. All those states are mixed states of positive and negative 
parities as in Eq. (f3l2|). We add a prime in the following to indicate each single-particle state 
has both positive-parity and negative-parity components. For example, a 0s[^ 2 state has 
S1/2 (positive-parity) and pi/ 2 (negative parity) components with an S1/2 component as a 
dominant one. For 14 the neutron orbits are filled up to 0p' 3 ^ 2 . The mixing probability 
of the opposite parity components for 0s' 1/2 is larger than 20% and that for 0p' 3/2 is a few 
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FIG. 8: (Color online) Kinetic energy (T) per particle and the sum of the potential energies 
from the central force Vc, the Coulomb force Vcoul> an d the three-body force V3B divided by mass 
numbers. The horizontal axis indicates mass numbers. The dashed lines correspond to the kinetic 
energy and the solid lines to the sum of the potential energies. The circle, triangle, and diamond 
symbols correspond to the Hartree-Fock (HF) calculation, the charge- and parity-projected Hartree- 
Fock (CPPHF) calculation with the normal tensor force (xt = 1-0), and the CPPHF calculation 
with the strong tensor force (xt = 1.5). 

% in 14 O. In the CPPHF method the tensor correlation energy is gained by parity mixing 
and, therefore, P m i x is a measure to indicate how each single-particle orbit contributes to 
the tensor correlation. The large P m i x for 0sy 2 indicates that this orbit is largely affected by 
the tensor correlation. In 16 O, the 0p^, 2 neutron orbit is added newly. The addition of the 
0p[/ 2 orbit reduces the parity mixing of the 0sy 2 because they have the same total angular 
momentum j = 1 / 2 . As the result, the tensor correlation energy becomes smaller. This effect 
is more significant for the strong tensor force case as seen in Fig. El In 22 the neutron 
0g?q/ 2 is filled. Because there are no j—5/2 below, P m ix's for the states filled already in 16 
do not change largely. In 24 the ls' x , 2 orbit is newly occupied. The occupation of ls'jy 2 
reduces P m i x 's for the states with j—1/2, 0s[, 2 and 0py 2 . The tensor correlation energy is 
enhanced, although the amount of the change is small. Finally, in 28 the neutron 0d' 3 , 2 is 
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FIG. 9: (Color online) Potential energies from the tensor force (Vt) and the LS force (Vls) divided 
by mass numbers. The horizontal axis indicates mass numbers. The solid lines correspond to the 
potential energy from the tensor force. The dashed line corresponds to the potential energy from 
the LS force in the CPPHF calculation with the normal tensor force. Other results for the LS force 
are almost unchanged. The meanings of the symbols are the same as in Fig. |HJ 

filled. Pmix's for the previously filled orbits change by this addition but the changes are not 
so large. 

The change of P m j x as shown above indicates that the large change of the tensor correlation 
energy from 14 to 16 in the strong tensor force case is caused by the blocking effect for 
the j=l/2 orbits. The blocking effect of this kind is shown to be important in the single- 
particle Zs-splitting in 5 He p]. The effect of the blocking is less significant for the excess 
neutron orbits. The effect of the blocking on binding energies in neutron excess oxygen 
isotopes seems to be small but the blocking may affect single-particle natures or collectivity 
in neutron excess oxygen isotopes, because the mixing probability is affected by the blocking 
effect. 

The negligible P m i x for the 0d' 5 ^ 2 orbit indicates that this orbit does not contribute to the 
tensor correlation although there is no other occupied orbits which have j=5/2. The main 
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FIG. 10: Probabilities of the mixing of the opposite parity components for single-particle states 
with neutron components as dominant ones in the charge- and parity projected Hartree-Fock 
calculation with the normal tensor force. The horizontal axis indicates mass numbers. The square, 
circle, triangle, asterisk, diamond, and star symbols correspond to the first j=l/2 (0s^ 2 ! s i/2 
dominant), the first j=3/2 (0p' 3 / 2 ', P3/2 dominant), the second j = 1 / 2 (0py 2 ; p 1 / 2 dominant), the 
first j=5/2 (0d' 5 ^ 2 ; cZ 5 y 2 dominant), the third j = 1 / 2 (lsy 2 ; s ly / 2 dominant), and the second j=3/2 
(Qdy 2 ; d 3 /2 dominant), respectively. 

part of the tensor correlation comes from the T = channel. Because the proton j—5/2 
orbit is not filled in neutron-excess oxygen isotopes, the 0d' 5 ^ 2 orbit is hard to contribute to 
the tensor correlation and the mixing probability of it becomes small. 

IV. SUMMARY 

We have studied the effect of the tensor force in the sub-closed-shell oxygen isotopes 
using the charge- and parity-projected Hartree-Fock (CPPHF) method. We have extended 
the CPPHF method to the cases with Hamiltonians including three-body forces, although 
the extension is straightforward. In the CPPHF method the parity and charge-number 
projections are performed before variation. By applying the CPPHF method to the oxygen 



22 



isotopes actually, we have found that a sizable potential energy from the tensor force is 
obtained in the CPPHF method while in the Hartree-Fock calculation quite a small potential 
energy from the tensor force is obtained. 

We have investigated 16 O in some details. The correlation energy from the tensor force is 
about 10 MeV for the normal tensor force case and about 40 MeV for the strong tensor force 
case. In the strong tensor force case the strength of the T\ ■ T2 channel in the tensor force is 
multiplied by 1.5. The opposite parity components induced in single-particle states by the 
tensor force have compact sizes as compared to the normal parity components. This indicates 
that for the tensor correlation high- momentum components are important as already found 

nnnnn 

in the alpha particle case p, [12J, [13J, [15J, [16|| . The present result infers that the importance of 



the high-momentum component for the tensor correlation is valid in heavier mass nuclei. We 
have also shown the density and the charge form factor calculated in the CPPHF method. 
The density in the CPPHF method is reduced around the center and is pulled in to the 
inside region. This is caused by the parity mixing and the shrinkage of single-particle wave 
functions of opposite parities. The effect of the shrinkage appears in the charge form factor as 
a tail in a high-momentum region, because the shrinkage of the single-particle wave functions 
induces high-momentum components in the density. 

In the results for the oxygen isotopes, the general tendencies for the binding energies 
and the matter radii are reproduced in the CPPHF calculation with the effective interaction 
adopted here, while the agreement of the binding energy with the experimental data is not 
so good compared to the Hartree-Fock-type calculations with recent effective interactions. 
The root-mean-square matter radii are well reproduced within error bars except for 14 O with 
both the normal tensor force and the strong tensor force. In the strong tensor force case, 
the matter radius of 14 becomes smaller and close to the experimental data. The reduction 
of the matter radius in 14 O with the strong tensor force is due to the shrinkage of single- 
particle wave function by the strong tensor correlation, which is large in 14 O. Because the 
Opi/2 neutron orbit is not occupied in 14 0, there are no blocking states for the Osi/2 proton 
orbit in the tensor correlation and, therefore, the tensor correlation becomes large. Actually, 
the correlation energy from the tensor force per particle amounts to more than 5 MeV in 
this case. For all the oxygen isotopes the calculated potential energy from the tensor force 
is in the same order as that from the LS force for the normal tensor force case. In the strong 
tensor force case it becomes about two times larger. In contrast to the potential energy 
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from the LS force, the potential energy from the tensor force in an LS-closed-shell nucleus 
does not become close to zero in the CPPHF calculation. In the Hartree-Fock calculation it 
becomes negligibly small because both the total spin and the total orbital angular momenta 
is almost zero in an LS-closed-shell nucleus. The sudden decrease of the potential energy 
from the tensor force from 14 to 16 is attributed to the blocking effect of the j—1/2 orbits. 
The blocking effect is also seen in the neutron-excess oxygen isotopes but does not affect 
the binding energy largely. 

In the present study we have applied the CPPHF method to the ground states of the 
sub-closed-shell oxygen isotopes. The application to odd-mass nuclei and open-shell nuclei 
to study the effect of the tensor force on single-particle natures and the change of collectivity 
by the tensor correlation in a neutron excess region are interesting because the tensor force 
changes the spin, the orbital angular momenta and the isospin of nucleon orbits simultane- 
ously, which is realized in the CPPHF method by the parity and charge mixing. As for an 
effective interaction, we combine the available effective interaction and the tensor force in 
the free space with some modifications. We need to use effective interactions which have the 
connection with the realistic nuclear forces to reveal the relation between nuclear structure 
and the underlying nuclear force and have the consistency between the tensor force and 
other forces like the central and LS forces. The study in such a direction is also important 
and now under progress. 
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